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Abstract — It is well-known that the behavior of many dy- 
namical processes running on networks is intimately related 
to the eigenvalue spectrum of the network. In this paper, we 
address the problem of inferring global information regarding 
the eigenvalue spectrum of a network from a set of local samples 
of its structure. In particular, we find explicit relationships 
between the so-called spectral moments of a graph and the 
presence of certain small subgraphs, also called motifs, in the 
network. Since the eigenvalues of the network have a direct 
influence on the network dynamical behavior, our result builds 
a bridge between local network measurements (i.e., the presence 
of small subgraphs) and global dynamical behavior (via the 
spectral moments). Furthermore, based on our result, we 
propose a novel decentralized scheme to compute the spectral 
moments of a network by aggregating local measurements of 
the network topology. Our final objective is to understand 
the relationships between the behavior of dynamical processes 
taking place in a large-scale complex network and its local 
topological properties. 

I. Introduction 

Research in complex networks has important applications 
in today's massive networked systems, including the Internet, 
the World-Wide Web (WWW), as well as social, biological 
and chemical networks [1]. The availability of massive 
databases, and reliable tools for data analysis provides a 
powerful framework to explore structural properties of large- 
scale networks [2]. In many real-world cases, it is impossible 
to efficiently retrieve and/or store the exact structure of a 
complex network due to, for example, a prohibitively large 
network size or privacy /security concerns. On the other hand, 
it is often possible to gather a great deal of information 
by examining local samples of the graph topology, such as 
the degree distribution, which are usually easy to collect. In 
this context, a challenging problem is to find relationships 
between the set of local samples of the network structure 
and its global functionality. 

Many dynamical processes on networks, such as random 
walks [3], virus/rumor spreading [4], [5], and synchro- 
nization of oscillators [6], [7], are interesting to study in 
the context of large-scale complex networks. The behavior 
of many of these processes are intimately related to the 
eigenvalue spectra of the underlying graph structure [8], [9]. 
For example, the speed of spreading of a virus is directly 
related with the spectral radius of the adjacency matrix of the 
network [4]. Hence, spectral graph theory provides us with 
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a framework to study the relationship between the network 
structure and its dynamical behavior 

Network motifs are small subgraphs that are present in a 
network with a much higher frequently than in random net- 
works with the same degree sequence [10]. There is both em- 
pirical and theoretical evidence showing that these subgraphs 
play a key role in the network's function and organization 
[11]. One of the main objectives of this paper is to explicitly 
relate global properties of a given network with the presence 
of certain small subgraphs that can be counted via local 
measurements. We focus our attention on global properties 
related with the network's eigenvalues, in particular, the so- 
called spectral moments. Since the spectral properties are 
known to have a direct influence on the network dynamical 
behavior, our result builds a bridge between local network 
measurements (i.e., the presence of small subgraphs) and 
global dynamical behavior (via the spectral moments). It 
is also worth pointing out that, although a set of spectral 
moments is not enough to completely describe the spectral 
distribution of a network, it allows us to extract a great deal 
of information. For example, Popescu and Bertsimas provide 
in [12] an optimization framework for computing optimal 
bounds on the properties of a distribution from moments 
constrains. More generally, there is a variety of techniques 
that can be applied to extract spectral information from a 
truncated sequence of spectral moments [8]. 

Based on our results, we also propose a novel decentral- 
ized algorithm to efficiently aggregate a set of local network 
measurements into global spectral moments. Our work is re- 
lated to [13], where a fully distributed algorithm is proposed 
to compute the full set of eigenvalues and eigenvectors of 
a matrix representing the network topology. In contrast, our 
approach is computationally much cheaper, since it does not 
require a complete eigenvalue decomposition. Furthermore, 
our approach also provides a clearer view of the role of 
certain subgraphs in the network's dynamical behavior 

The rest of this paper is organized as follows. In the 
next section, we review graph-theoretical terminology and 
introduce definitions needed in our derivations. In SectionHID 
we derive explicit relationships between the moments of the 
eigenvalue spectrum and local network measurements. Based 
on these expressions, we introduce a distributed algorithm 
to compute these moments in Section |IV] We conclude the 
paper mentioning some future work. 

II. Definitions & Notation 

Let Q = (V, £) denote a simple undirected graph (with no 
self-loops) on n nodes, where V {Q) = {vi, . . . , u„} denotes 



the set of nodes and £ (G) C V (G) x V (G) is the set of 
undirected edges. If {vi,Vj} e £ (G) we call nodes Vi and 
Vj adjacent (or neighbors), which we denote by Vi ^ Vj. 
The set of all nodes adjacent to a node u G V (G) constitutes 
the neighborhood of node v, defined by TV" = {w E V (G) ■ 
{v, w} G £ (G)}, and the number of those neighbors is called 
the degree of node v, denoted by degv or d„. 

We define a walk of length k from vq to Vk to be an ordered 
sequence of nodes {vo,vi, ...,Vk) such that Vi ^ Vi+i for 
i = 0,1, ...,fc — 1. If Vq = Vk, then the walk is closed. 
A closed walk with no repeated nodes (with the exception 
of the first and last nodes) is called a cycle. Triangles and 
quadrangles are cycles of length three and four, respectively. 
We say that a graph G is connected if there exists a walk 
between every pair of nodes. Let d (u, w) denote the distance 
between two nodes v and w, i.e., the minimum length of 
a walk from v to w. We define the diameter of a graph, 
denoted by diam {G), as the maximum distance between 
any pair of nodes in G- We say that v and w are fc-th 
order neighbors if d (u, w) = fc, and define the fc-th order 
neighborhood of a node v as the set of nodes within a 
distance fc from v, i.e., Af^ = {w G V {G) ■ d {v, w) < k}. A 
fc-th order neighborhood, induces a subgraph G^ ^ G with 
node-set Af^ and edge-set £J! defined as the subset of edges 
of £ (G) that connect two nodes in A/'^y. 

Graphs can be algebraically represented via the adjacency 
matrix. The adjacency matrix of an undirected graph G, 
denoted by Ag = [a^j], is an n x n symmetric matrix defined 
entry-wise as = 1 if nodes Vi and Vj are adjacent, and 
Qij = otherwis^ll 

III. Spectral Analysis via Subgraph Embedding 

In this section, we derive an explicit relationship between 
the spectral moments of the adjacency matrix and the pres- 
ence of certain subgraphs in G- We say that a graph "H is 
embedded in if H is isomorphicjl to a subgraph in G- 
The embedding frequency of "H in G, denoted by F{T-L, G), 
is the number of different subgraphs of G to which H is 
isomorphic. The term network motif is used to designate 
those subgraphs of G that occur with embedding frequencies 
far higher than in random networks with the same degree 
sequence [10]. Theoretical and experimental evidence shows 
that some of these motifs carry significant information about 
the network's function and organization [11]. In this section, 
we derive an explicit expression for the spectral moments 
as a linear combination of the embedding frequencies of 
certain subgraphs. Our results provide a direct relationship 
between the presence of network motifs and global properties 
of the network, in particular, spectral moments. In the coming 
subsections, we first provide a theoretical foundation for 
our analysis. Second, we derive explicit expressions for the 
spectral moments in terms of network metrics, such as the 

'For simple graphs with no self-loops, an = for all i. 

^An isomorphism of graphs Q and is a bijection between the vertex 
sets V{G) and V(W), / : V{G) -+ V('H), such that any two vertices u and 
V of Q are adjacenct in G if and only if f{u) and f{v) are adjacent in H. 



degree sequence or the number of triangles in the graph. 
In the next Section, we propose a decentralized algorithm 
to compute the spectral moments of a network based on 
decentraUzed subgraph counting. 

A. Spectral Moments and Subgraph Embedding 

Algebraic graph theory provide us with the tools to relate 
topological properties of a graph with its spectral properties. 
In particular, we are interested in studying the spectral 
moments of the adjacency matrix of a given graph. We 
denote by the set of eigenvalues of Ag and define 

the fc-th spectral moment of the adjacency matrix as 



rrik (Ag) 



1 " 



(1) 



From algebraic graph theory, we have the following result 
relating the fc-th spectral moment of a graph G with the 
number of closed walks of length k in G [14]: 

Lemma 1: Let be a simple graph. The fc-th spectral 
moment of the adjacency matrix of G can be written as 



mk{Ag) = - 



(fc) 



(2) 



is) 

where g is the set of all possible closed walks of length 
fc in dl 

Corollary 2: Let be a simple graph. Denote by Eg and 
Ag the number of edges and triangles in G, respectively. 
Then, 

mi{Ag) ^ 0, m2{Ag) = 2Eg/n, and m3{Ag) = 6 Ag/n. 

(3) 

In the following, we develop on the above Lemma to 
derive an expression for nik {Ag) , for any fc, as a linear 
combination of the embedding frequencies of certain sub- 
graphs. Although the first 3 moments present very simple 
expressions in Corollary [2] the larger the fc, the more 
involved those expressions become. In what follows, we 
describe a systematic procedure to derive these expressions 
efficiently. 

First, we need to introduce some nomenclature. Given 
a walk of length k in G, w ^ (wq, wi, Wfc), we denote 
its node-set as V (w) — {vQ,...,Vk} and its edge-set as 
£ (w) = U^^i {vi^i,Vi}. Hence, we define the underlying 
simple graph of a walk w as H (w) = (V (w) , £ (w)). We 
say that a simple subgraph h C G spans the walk w if ft, is 
isomorphic to H (w). Notice how different walks can share 
the same underlying simple graph. We also denote by I (h) 
the unlabeled simple graph isomorphic to a given graph h. 
Applying the function 2 to the underlying graph of a given 
walk w, we obtain an unlabeled graph g ~ I {H {w)). Notice 
how different walks, not even sharing the same edge-set, 
can share the same unlabeled underlying simple graph. For 
example, any closed triangular walk w-s ~ {vi,Vj,Vk,Vi) 

'We denote by \Z\ the cardinality of a set Z. 



with Vi 



Vi gives rise to the same unlabeled 



graph I {H (103)), namely, an unlabeled triangle. 

We denote the set of closed walks of length k in Q 
sharing the same underlying simple graph h as TZh = 



{ 



w £^ 



P s.t. n {w) 



h 



|. Inversely, consider the set of 



>'tv><^o 



of 



Fig. 1 . In this table, we represent tlie set of nonisomorpliic connected grapli 
witli at most 4 nodes and 4 edges. The set of subgraphs in I4 are those 



the computation of the 4-th spectral moment. 



Hence, we can rewrite the above as 



all closed walks of length k in Q, '^g . We can define 

another set, which we denote by H^, containing all pos- with coiresponding coefficient c.^) > 0. These subgraphs are" involved in 

sible simple graphs spanning walks in "^g , i.e., = 

{wk) ,Wk G ^'tf'!' Furthermore, we define the set of 
unlabeled graph associated with walks of length k as Ij, — 

{71 (wk)) ,Wk G In what follows, we derive an 

expression for the fc-th spectral moment as a linear combina- 
tion of the embedding frequencies of the (unlabeled) graphs 
in Ifc. 

Note that the mapping H : Hk, is a surjection 

onto Hfe that induces a partition in ^^g \ namely, each graph 



(fe) 



geifc heSg 



/i G Hfc (with labeled nodes), induces a partition subset ^^^'"^ ^^^^ ^^^'^ the fact that \Rh 



,(fc) 



in the last 



(k) 

T^h Q ■ Similarly, the mapping I : — > I^- is 
a surjection onto 1^ and induces a partition in H^-. Each 
unlabeled graph g G Ife defines a partition subset Sg C 
defined as Sg ~ {/i G Hfe s.t. I (h) = g}. Note that for 
g G Ik, \Sg \ is the number of subgraphs of Q isomorphic to g; 
hence, \Sg\ equals the embedding frequency F{g,Q). This 
observation shall be important in the next Section, where 
we design a decentralized algorithm to compute spectral 
moments. 

Based on the above, we derive the following result regard- 
ing the spectral moments of Q: 

Theorem 3: Given a simple graph Q, its fc-th spectral 
moment can be written as 
1 



equality. Therefore, substituting the above in (|2]i we obtain 
the expression for the spectral moments in the statement of 
the theorem. ■ 
The expression in the right-hand side of (|4|i is a linear 
combination of the embedding frequencies of the set of 
subgraphs g G Ifc. The coefficients uJk (g) in this Unear 
combination are defined as the number of closed walks of 
length k with unlabeled underlying graph g. In what follows, 
we describe an algorithm to determine the set Ife and compute 



,(fe) 



In order to compute the first K spectral moments, we 



mk{AQ 



■E 



F{g.Q). 



where uj. 



(fe) 

9 



[we^f^ ■.T{n{w))^g] 



(4) 



in other 



words, ujg is the number of all possible closed walks of 
length k m g with underlying unlabeled graphs isomorphic 
to g. 

Proof: Consider the composition function J = I o 



H. Since J : * 



(k) 



(k) 

have to study the sets of closed walks "i!g ' , for k < K. For 
each particular value of fc, the set Ife are unlabeled connected 
graphs with at most k edges. Also, the maximum number of 
nodes visited by walks in '^'g'^ are equal to k. We denote 
by Gfe the set of all (unlabeled) nonisomorphic connected 
graphs with at most k nodes and at most k edge^; hence, 
Ifc C Gfe. For example, we illustrate all the graphs in G4 
(excepting the isolated-node graph) in Fig. 1. 

In order to compute the fc-th moment via (HI, we need to 
compute CLig*^' for all g G Ife, where uj^f^^ can be interpreted as 
the the number of closed walks of length k in g that visit all 

the nodes and edges of q (hence, its underlying simple graph 

(fc) 



Ifc is a composition of surjective is .9)- The computational complexity of computing is the 



functions, it is itself a surjective function. Thus, J induces 
a trivial partition in ^I'g'^' as follows. Given an unlabeled 
graph g G Ife, we define a partition subset Tg C ^I'^'"'' as 
Tg = [we ^^g'^ s.t. J{w) = .9} = Uftes, Rh- Based on 
this partition, we can compute the number of closed walks 
of length k in Q, Wg \ as 



E i^.i = E E 1^'' 



Note that, for a particular (labeled) graph h G Hfc, the 
number of closed walks of length k such that their underlying 
simple graph is h, i.e., is a quantity that depends 

exclusively on the topology of the unlabeled version of h, 
i.e., g = We define this quantity as ujk{g) = \Rh\- 



same as the one of counting the number of Eulerian walks in 
an undirected (multi)graph. This counting problem is known 

to be #P-complete [15]. Hence, there is not a closed-form 

(fc) 

expression for ujg , and it has to be computed via a brute- 
force combinatorial exploration over all the possible closed 
walks of length k visiting all the nodes and edges of the 
subgraph g. This exploration can be performed in reasonable 
time for subgraphs g of small and medium size - which 
are the ones we are interested in. Although computationally 

expensive, this computation is done once and for all for 

(fc) 

a particular g. In other words, once ujg is computed, the 
coefficient for the spectral moment in ^ are known for any 
given Q. 

''This set can be easily determined using the command ListGraphs 
included in the package Combinatorica included in Mathematica. 
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Fig. 2. In the above table, we represent the set of nonisomorphic connected graphs with at most 7 nodes and 7 edges. The set of subgraphs in I5, Ig 
and I7 are those with corresponding coefficients ujg^^ , ojg*' , and uig"^^ > 0, respectively. 



(k) 

For convenience, we provide the coefficients uig for all 4-cycles as: 
the graphs in G4 in Fig. 1. Using this table as an example, 

we also observe that those graphs t G such that t ^ 1^, m4{Ag) 
have wp' = (since cjj'''' being zero indicates that there is 
no walk in t satisfying the conditions to be part of the set 
Ife). For convenience, we also provide a list of the graphs 



1 r 



(4) 



2W2 -Wi+ 8$^ 



(6) 



(k) 

in Ua:=,5^6.7Ia; and the associated coefficients uig in Fig. 2 
(where we have left blank those cells where ujg = 0). In the 
following paragraphs, we illustrate the usage of Theorem |3] 
to derive explicit expressions for the first 5 spectral moments 
of a given graph Q. 

From the subgraphs and coefficients in Fig. 1, we can 
compute the 4*'' spectral moment via ^ as follows 



m4(Ag) = - 

n 



1 r 



2Eg + AAg 



,(4) 



(5) 



where Eg is the number of edges, Ag is the number of 
wedge-graph^ in Q, and ^'g ^ is the number of 4-cycles in 
Q (see subgraphs and coefficients in Fig. 1). Thus, in order to 
compute the 4*'' spectral moment, we must be able to count 
the number of wedge-graphs and 4-cycles in Q. Although 
the number of edges and cycles in a graph are common 
network metrics, the number of wedges is not. It is then 
convenient to rewrite the number of wedges in terms of more 
familiar network metrics. In fact, we can rewrite the number 
of wedges in terms of the degree sequence of the graph as 
follows; 



degj 
2 



where Wr = J2v=i i^^S "vf is the r-th power-sums of the 
degree sequence {degw}"^^. Since Eg = |Z^"=idcgu in 
a simple graph, we can write |5] in terms of power-sums and 

wedge graph is isomorphic to a chain graph of length 2. 



We illustrate this result in the following example. 

Example 4: Consider the n-ring graph i?„ (without self- 
loops). We know that the eigenvalues of the adjacency 
matrix of the ring graph A/j^ are 1 2 cos i ^ } ■ Hence, 
the 4-th moment is given by the summation 7714 (Aij^) = 
i Y^^Zo (2cosi^) , which (after some computations) can 
be found to be equal to 6 for n ^ {2,4}. We can apply 
(|5]i to easily reach the same result without performing an 
eigenvalue decomposition, as follows. In the ring graph, we 
have Eg = n, A^,^ = n for n 7^ 2, A^^^ = 0, = for 
n ^ 4, and = 1. Hence, we obtain the same value for 
mi{AR^) directly from (|5]l. 

From the subgraphs in Fig. 2, and the row of coefficients 



,(5) 



we have the following expression for the 5*'' spectral 



moment via dUi; 

m^{Ag) 



30Ag + lOTg + 10$ 



(5) 



(7) 



where Ag and is the number of triangles and 5-cycles 
in Q. Also, we represent by Tg the number of subgraphs 
isomorphic to the sixth subgraph in the top row of Fig. 
2, counting from the left. Analyzing the structure of this 
subgraph, we have that Tg = Y^l^iidi — 2)ls!g\ where 
Ag'' is the number of triangles in Q touching node i. 
Hence, defining the clustering-degree correlation coefficient 
as C\d ~ Sr=i ^j^e'' we can write Tg ~ CacI — 6Ag, and 
the 5*'' spectral moment as: 



m5{Ag) 



1 



10$^^' - 30Ag + lOCAd 



(8) 



Notice that, apart from triangles and 5-cycles, the clustering- 
degree correlation influences the 5-th moment; hence, non- 
trivial variations of the local clustering with the degree, as 
reported in [16], are relevant in the computation of the 5-th 
spectral moments. 

Remark 5: The main advantage of (|4]i in Theorem [3] may 
not be apparent in graphs with simple, regular structure. 
For these graphs, an explicit eigenvalue decomposition is 
available and there may be no need to look for alternative 
ways to compute spectral moments. On the other hand, in 
the case of large-scale complex networks, the structure of 
the network is usually very intricate, in many cases not even 
known exactly, and an explicit eigenvalue decomposition can 
be very challenging, if not impossible. It is in these cases 
where the alternative approach proposed in this paper is most 
useful. 

As we show in the next section, the spectral moments 
can be efficiently computed via a distributed approach from 
aggregation of local samples of the graph topology, without 
knowing the complete structure of the network. 



IV. Distributed Computation of Spectral 
Moments 

According to (|4]i, the fc-th spectral moment is equal to 
a linear combination of the embedding frequencies F {g,Q) 
for g e Ife. On the other hand, computing the embedding 
frequencies in large-scale networks can be challenging in 
a centralized framework, since the computational cost of 
counting subgraphs rapidly grows with the network size. In 
this section, we introduce an efficient decentralized approach 
to compute the embedding frequencies of subgraphs from 
local samples of the network topology. 

During this section, we assume that there is an agent 
in each node v G V (G) that is able to access its r- 
th neighborhood, Q^!. A naive approach to compute the 
embedding frequency of a particular subgraph g would be 
to compute the embedding frequency of g in each neighbor- 
hood, F{g,Q^), and sum them up. This approach obviously 
does not work because this particular subgraph g Q Q can 
be in the intersection of multiple neighborhoods; hence, that 
subgraph would be counted multiple times. In what follows, 
we propose a decentralized counting procedure that allow us 
to know how many times a particular subgraph is counted. 

First, we need to introduce several definitions: 

Definition 6: We say that a graphical property Pg is 
locally measurable within a radius r around a node v if 
Pg is a function of A/";, i.e., Pg ^ f {J\f^). 

Definition 7: We say that a subgraph h C Q is locally 
countable within a radius r around a node v if h C_ . 

For example, both edges and triangles touching a node v 
are locally countable within a radius 1. Also, the number of 
quadrangles touching v is locally countable within a radius 
2. Furthermore, a wedge is locally countable within a radius 
1 only by the node at the center of the wedge, but it is not 
countable by the nodes at the extremes of the wedge. On the 
other hand, the wedge is locally countable within a radius 2 



by all the nodes in the wedge. In this examples, we observe 
how not all the nodes being part of a subgraph h C Q have 
to be able to locally count the subgraph. In particular, if the 
radius r is smaller than \diam (h) /2], none of the nodes in 
h is able to count the subgraph locally. On the contrary, if the 
radius r is greater or equal to the diameter of the subgraph, 
all the nodes in h are able to locally count it. In the middle 
range, \diam (h) /2] < r < diam (h), some nodes are able 
to locally count h and some are not. 

Definition 8: For a given value of r, the set of detector 
nodes of a given subgraph h C Q is defined as 

z?!;^ = {v ev{h) s.t. hcg:;:}. 

Note that, although there can be other nodes u ^ V (h) 
able to locally count h, we do not include them in the set of 
detector nodes Dh . Also note that, given an unlabeled graph 

(r) 

g, the size of Df^/ for all the subgraphs hi C Q isomorphic to 
g depends exclusively on the structure of g. In other words. 



we have that 



D 



given g, the value of 
as follows: 



(r) 



D 



(r) 



for all i — I, F {g, Q). For a 
can be algorithmically computed 



1) Initialize M := 0; 

2) For each node u G V (g); 



3) 
4) 
5) 
6) 



Compute 5g{u) := maxy^g d{u,v); 
If Sg{u) < r, Tlien M := M + 1; 



End For, 



D 



(r) 



M. 



D 



(r) 



for g e 



For convenience, we include the values of 
Ufe=4.5,6l/£ for radius ?■ = 1,2, and 3 in Fig. 4. 

After these preliminary results, we now describe a novel 
algorithm to distributedly compute F (g, Q). Let us consider 
a particular node w G V {Q). We define the local embedding 
frequency of subgraph g within a neighborhood of radius 
r around node v, denoted by H {g,J\f^), as the number of 
different subgraphs li C isomorphic to g such that v € 
V{li). Our distributed algorithm is based on the following 
result: 

Theorem 9: Let Q he a simple graph. The spectral mo- 
ments of the adjacency matrix of Q can be written as 



rukiAg) 



n ^ — ' 



E 



wfe (g) 



D 



(r) 



H{g,Kn 



(9) 



vev{g) geifc l^g 
Proof: Note that the local embedding frequency only 
count subgraphs li touching node v; hence, H{g,Af^) < 
F{g,J\f^). Furthermore, we have that 



E 



H{g,N:)^\Dg\ F{g,g) 



(10) 



13, 



since every graph h C Q isomorphic to g is counted 
times in the above summation. Therefore, substituting ( [TOl 
in dill, we obtain the statement of the Theorem. ■ 

Based on Theorem |9] it is straightforward to compute the 
fc-th spectral moment distributedly as follows. First, define 







r=l 



r=3 



Fig. 3. In the above table, we represent several nonisomorphic connected graphs, and the con'esponsing values of \Dg for r G {1, 2, 2} 



the following local variables 



E 

geik 



wfc (g) 



(r) 



for all f G V{Q). Note that /xjj'' (u) is locally measurable 
within a radius r around each node v, namely, each agent 
is able to distributedly compute /iJJ (v) for all v G V{Q). 
Thus, from (|9]l, we have that the fc-th spectral moment can 
be computed via a simple distributed averaging of /ijl (v). 

Remark 10: The maximum order of the spectral moment 
that can be computed via this distributed approach depends 
on the radius r of the neighborhood that each agent can 
reach. In particular, in order to compute the fc-th moment, 
we should be able to detect all the graphs g G Ifc- It is easy 
to prove that the fc-ring with diameter [fc/2j is the graph in 
Ifc with the maximum diameter. As we mentioned before, in 
order for a particular subgraph h to be locally countable, the 
radius r must be greater or equal to \diara (h) /2]. Hence, 
for a particular r, we can distributedly compute moments up 
to an order fc,nax = 2r + 1. 

V. Conclusions and Future Research 

In this paper, we have derived explicit relationships be- 
tween spectral properties of a network and the presence of 
certain subgraphs. In particular, we are able to express the 
spectral moments as a linear combination of the embedding 
frequencies of these subgraphs. Since the spectral properties 
are known to have a direct influence on the dynamics 
of the network, our result builds a bridge between local 
network measurements (i.e., the presence of small subgraphs) 
and global dynamical behavior (via the spectral moments). 
Furthermore, based on our result, we have also developed a 
novel decentralized algorithm to efficiently aggregate a set of 
local network measurements into global spectral moments. 
Our approach is based on a an efficient decentralized ap- 
proach to compute the embedding frequencies of subgraphs 
from local samples of the network topology. 

Future work involves to extend our methodology to di- 
rected graphs and graphs with self-loops, such as those ap- 
pearing in transcription networks. Also, we are interested in 
developing techniques to extract explicit information regard- 
ing the dynamical behavior of a network from a sequence 



of spectral moments. Furthermore, it would be interesting 
to find a fully decentralized algorithm to iteratively modify 
the structure of a network of agents in order to control its 
dynamical behavior. We find particularly interesting the case 
in which individual agents have knowledge of their local 
network structure only (i.e., myopic information), while they 
try to collectively aggregate these local pieces of information 
to find the most beneficial modification of the network 
structure. 
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